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ABSTRACT 

We combine new high resolution imaging and spectroscopy from Keck/NIRC2, Discovery Chan¬ 
nel Telescope/DSSI, and Keck/HIRES with published astrometry and radial velocities to measure 
individual masses and orbital elements of the GJ 3305 AB system, a young (^20 Myr) M-l-M binary 
(unresolved spectral type MO) member of the /3 Pictoris moving group comoving with the imaged 
exoplanet host 51Eri. We measure a total system mass of 1.11 ± 0.04 Mq, a period of 29.03 ± 0.50 
year, a semimajor axis of 9.78 ± 0.14 AU, and an eccentricity of 0.19 ± 0.02. The primary component 
has a dynamical mass of 0.67 ± 0.05 Mq and the secondary has a mass of 0.44 ± 0.05 Mq. The 
recently updated BHAC15 models are consistent with the masses of both stars to within 1.5cr. Given 
the observed masses the models predict an age of the GJ 3305 AB system of 37 ± 9 Myr. Based on the 
the observed system architecture and our dynamical mass measurement, it is unlikely that the orbit 
of 51Erib has been significantly altered by the Kozai-Lidov mechanism. 

Subject headings: astrometry — binaries: close — stars: fundamental parameters — stars: individual 
(GJ3305AB) 


1. INTRODUCTION 

Loose associations of young, nearby (<70 pc) stars 
with common ages, kinematics, and origins have been 
a subject of increasing interest (Zuckerman et al. 2004; 
Shkolnik et al. 2012; Malo et al. 2013). Because of their 
proximity to Earth, these young moving groups (YMGs) 
are excellent targets to study pre-main sequence (PMS) 
stellar and substellar evolution, protoplanetary and de¬ 
bris disk structure, and giant planet formation at ages 
between distant star-forming regions and old field stars 
(e.g. Close et al. 2005; Nielsen & Close 2010). About 10 
YMGs containing hundreds of objects between 8 and 120 
million years old are known (e.g. Torres et al. 2008). 

As these moving groups are amenable to numerous age 
dating methods, including kinematic techniques, they 
provide the opportunity to measure dynamical masses 
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of PMS low-mass binary objects and test stellar evolu¬ 
tion models (Stassun et al. 2014). Generally, PMS stel¬ 
lar masses are inferred by comparing a star’s temper¬ 
ature, luminosity and metallicity to model predictions 
(e.g. Schaefer et al. 2014). These models are poorly 
constrained by observations and may induce systematic 
offsets (Dupuy et al. 2009, 2014). Worse yet, different 
models predict disparate masses, primarily due to un¬ 
certainties in the treatment of convection in low-gravity 
atmospheres (Baraffe et al. 2002), stellar accretion his¬ 
tory (Baraffe & Ghabrier 2010), and molecular line lists 
(Baraffe et al. 2015). In some cases, model-predicted 
masses can differ by a factor of two or more (Hillenbrand 
& White 2004; Schlieder et al. 2014). Dynamical mass 
measurements of binary stars with known ages are essen¬ 
tial to test models. 

Recently, Macintosh et al. (2015) presented 51Erib, 
the first exoplanet discovery from the Gemini Planet Im¬ 
ager. The planet has a mass of « 2 Mjup (assuming a hot 
start model), a projected separation of 13 AU, a temper¬ 
ature of 600-750 K, and a T4.5-T6 spectral type. GJ 3305 
is known to be a binary with combined spectral type MO 
(Kasper et al. 2007). Feigelson et al. (2006) identified 
GJ 3305 and 51 Eri as an FO-MO common proper motion 
pair, separated by 66" or ~2000 AU. 

As a binary system, a dynamical mass can be measured 
for both stars in GJ 3305 AB. As both stars are members 
of the /? Pictoris moving group, an approximate age of 
the system is known (24 ± 3 Myr; Binks & Jeffries 2014; 
Mamajek & Bell 2014; Bell et al. 2015). While most 
dynamical masses of M dwarfs are limited by distance 
uncertainties, 51 Eri has a parallax from Hipparcos mea¬ 
sured to a precision of 1%. Combining this parallax with 
15 years of imaging and RV data enables us to determine 
the system orbital parameters, elucidating the architec¬ 
ture of this 4—or more—body system. 

In this paper, we combine RV and astrometric observa- 


2 


tions of GJ 3305 AB to measure orbital parameters and 
masses for each star. We compare these masses to model 
predictions and discuss the possible implications of this 
binary pair on the long-term evolution of the orbit of 
51 Erib. 


2. DATA COLLECTION AND REDUCTION 

GJ 3305 AB has been imaged and resolved many times 
(Kasper et al. 2007; Bergfors et al. 2010; Delorme et al. 
2012; Janson et al. 2012, 2014). The system was also im¬ 
aged with NIRG2 (Wizinowich et al. 2000) in one unpub¬ 
lished epoch in 2001 available in the Keck Observatory 
Archives (KOA, PI Zuckerman). In this work, we com¬ 
bine these data with hve observations from 2002 to 2015, 
three using Keck/NIRC2 and one with the Differential 
Speckle Survey Instrument (DSSI, Horch et al. 2009) at 
the Discovery Channel Telescope at Lowell Observatory. 

All NIRC2 data were obtained with the narrow cam¬ 
era mode, which has a field of view of 10'.'2 x I 0'.'2 and 
a plate scale of 9.952 mas pixel“^ (Yelda et al. 2010). 
All images were flat fielded and cleaned of bad pixels 
and cosmic rays. Astrometry and relative photometry of 
GJ 3305 was derived by simultaneously htting three bi¬ 
variate Gaussians to each component following Liu et al. 
( 2010 ). 

DSSI allows for simultaneous observations in two fil¬ 
ters. We use the DSSI R and I filters, with central wave¬ 
lengths 692 and 880 nm and FWHMs of 40 and 50 nm. 
We obtained 1000 40-ms exposures in each channel simul¬ 
taneously. The data were then reduced following Horch 
et al. (2015). Specifically, the autocorrelation of each 
frame was calculated and summed over all exposures, and 
the near-axis subplanes of the image bispectrum were 
calculated. To create a reconstructed image, the Fourier 
transform of the autocorrelation of the binary was di¬ 
vided by that of a nearby point source (HR 1415). The 
square root of this value is taken, and the result com¬ 
bined with a phase function derived from the bispectral 
subplanes. The pixel scale (19 mas pixel“^ in R and 
20 mas pixel”^ in I) and orientation of the detector were 
found by observing several widely separated binaries with 
known astrometry. Our astrometry is listed in Table 1. 

The GJ 3305 binary system has also been monitored 
spectroscopically. One Keck/HIRES spectrum from 2003 
exists in the KOA (PI Zuckerman); we measure the RV 
following Kraus et al. (2015). We combine this spectrum 
with nine additional spectra from Bailey et al. (2012), 
Shkolnik et al. (2012), and Elliott et al. (2014). In all 
cases, the RVs were calculated treating the system as an 
SBl. We take the reported RV and uncertainty for each 
observation, but assume the flux from the secondary is 
non-negligible, as explained in Section 3. 

3. ANALYSIS 

We infer the orbital parameters of GJ 3305 AB by com¬ 
paring the astrometric and RV data to a Keplerian orbit 
model at each of the observation times. A parallax, as¬ 
trometric orbit, and SBl RV data can be combined to 
measure individual masses of each star (e.g. Bean et al. 
2007). There is no measured parallax for GJ 3305, so 
we adopt the Hipparcos distance to 51 Eri/HIP 21547: 
29.43 ± 0.30 pc (van Leeuwen 2007). These two comov¬ 
ing systems have a projected separation of 1940 ± 20 
AU, or 0.01 pc. It is unlikely that the radial distance 


between the two could be significantly larger while re¬ 
maining bound; we apply this parallax as a prior on the 
distance to GJ 3305. 

We then fit for nine additional parameters that define 
the orbits of the two stars as viewed from Earth. Of 
these, seven can be obtained from astrometry. These 
parameters are the eccentricity vectors y/e cos uj and 
•Resina;, the time of periapse fp, the period P, the total 
mass Ml -|- M 2 , the inclination i, and the longitude of 
the ascending node SI. We parameterize the eccentricity 
vector in this manner following Eastman et al. (2013). 
The RV data can provide additional information about 
several of these (not Mi + M 2 or SJ directly), also allow¬ 
ing us to fit the systemic RV 7 and the secondary mass 
M2. 

We include ten additional terms to account for possi¬ 
ble systematics in the datas. This star has been imaged, 
resolved, and published by four different groups. We 
account for the possibility each group may have under¬ 
estimated their uncertainties on the orbital separation 
and position angle by a multiplicative factor by includ¬ 
ing a systematic error term on the measured positions 
from each group, allowing outlier points to be down¬ 
weighted without manually choosing specific points to 
downweight. We do the same with our reductions of 
both archival and new data, allowing for separate sys¬ 
tematic error terms on our data from Keck/NIRG2 and 
DGT/DSSI, providing a total of six systematic error 
terms. We allow the uncertainties on each dataset to 
be inflated up to a factor of five. 

Similarly, we allow for the possibility that the uncer¬ 
tainties in the RVs may be underestimated, possibly due 
to stellar variability (Moulds et al. 2013), errors in sys¬ 
temic RVs of standard stars, or drifts in the stability of 
the spectrographs. As our RV data originate from four 
sources, we allow each to have its own systematic error 
term, analogous to the jitter term commonly applied in 
RV orbit fits of exoplanets (e.g. Johnson et al. 2011): 


log£ (X 
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( 1 ) 

Here, C is the likelihood of the data given some under¬ 
lying physical model, CTo.z is the observed uncertainty on 
the zth data point, Us the systematic error associated 
with each particular set of observations, fi{t) the RV 
model evaluated at time t, and Vi{t) the observed RV at 
each t. Maximum likelihood jitter values range from 0.13 
km s“^ for the 2003 HIRES data to 0.57 km s“^ for the 
UVES data, suggesting stellar jitter is significant in the 
RV data, as expected for young stars. 

In all cases, one set of lines are observed because the 
RV separation is smaller than the line width. We ex¬ 
pect each RV measurement to be the flux-weighted sum 
of the two individual RVs. At each step, we calculate 
the RVs for each star, weighting them according to their 
expected flux contribution in each bandpass, using the 
observed flux ratios in the visible and near-IR as priors 
and assuming an additional 0.1 mag of variability in the 
optical and 0.05 mag in the near-IR. 

We neglect the possibility that 51 Eri could contribute 
signihcantly to the observed RV signal. Following Equa¬ 
tion 1 of Montet et al. (2015), the maximum RV acceler- 
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TABLE 1 

Data for GJ 3305 AB 


Epoch 

Bandpass 

RV 

Contrast 

Separation 

Position Angle 

Source 

(Year) 


(km s“^) 

(A mag) 

(mas) 

(deg) 


2001.910 

H2(^'=1-0) 


1.00 ±0.02 

286 ± 1 

198.1 ± 0.1 

This Work 

2002.162 

H 


1.02 ±0.02 

275.4 ± 1.5 

197.9 ±0.2 

This Work 

2003.05 

K 


0.94 ± 0.05 

225 ± 5 I 

195.0 ± I. 5 I 

Kasper et al. (2007) 

2003.195 

H 


0.99 ±0.01 

217 ± 1 

196.8 ±0.1 

This Work 

2004.02 

L' 



159 ±2 

194 ± 1 

Delorme et al. (2012) 

2004.95 

L' 


0.88 ±0.28 

93 ±2 

189.5 ± 0.4 

Kasper et al. (2007) 

2008.88 

SDSS 


1.39 ±0.16 

218 ±2 

20.3 ±0.3 

Bergfors et al. (2010) 

2008.88 

SDSS i' 


2.57 ± 0.05 

218 ±2 

20.3 ±0.3 

Bergfors et al. (2010) 

2009.13 

SDSS i'+z' 



231 ± 2 

19.2 ±0.3 

Janson et al. (2012) 

2009.90 

L' 



269 ±3 

18.6 ± 1.0 

Delorme et al. (2012) 

2009.98 

L' 



272 ±3 

19.2 ± 1.0 

Delorme et al. (2012) 

2010.10 

SDSS V 


1.34 ± 0.01 

284 ± 3 

18.5 ±0.6 

Janson et al. (2012) 

2010.10 

SDSS i' 


3.73 ±0.01 



Janson et al. (2012) 

2010.81 

SDSS V 



297 ± 3 

19.4 ± 0.3 

Janson et al. (2014) 

2011.67 

L' 



303 ±3 

18.1 ± 1.0 

Delorme et al. (2012) 

2011.87 

SDSS V 



295 ±4 

18.5 ±0.3 

Janson et al. (2014) 

2012.01 

SDSS V 



307 ± 3 

18.2 ±0.3 

Janson et al. (2014) 

2014.629 

Bry 


0.92 ±0.01 

244 ± 1 

16.8 ±0.1 

This Work 

2014.746 

DSSI R 


1.89 ±0.04 

239 ± 1 

16.4 ±0.2 

This Work 

2014.746 

DSSI I 


1.17 ± 0.03 

240 ± 1 

16.1 ± 0.2 

This Work 

2015.653 

K 


0.93 ±0.01 

199 ± 1 

15.6 ±0.1 

This Work 

2015.653 

H 


0.99 ±0.01 

198 ± 1 

15.6 ±0.1 

This Work 

2015.653 

J 


0.97 ± 0.01 

199 ± 1 

15.6 ±0.2 

This Work 

2015.653 

Y 


1.06 ±0.03 

200 ± 1 

15.6 ±0.1 

This Work 

2003.796 

HIRES V 

19.41 ± 0.38 




This work 

2004.884 

NIRSPEG K 

19.86 ±0.05 




Bailey et al. (2012) 

2005.862 

NIRSPEG K 

20.55 ±0.06 




Bailey et al. (2012) 

2005.971 

HIRES V 

21.70 ±0.30 




Shkolnik et al. (2012) 

2006.014 

NIRSPEG K 

20.82 ± 0.05 




Bailey et al. (2012) 

2006.016 

NIRSPEG K 

20.95 ±0.05 




Bailey et al. (2012) 

2006.019 

NIRSPEG K 

20.95 ± 0.05 




Bailey et al. (2012) 

2011.778 

UVES Blue 

24.40 ± 0.04 




Elliott et al. (2014) 

2001.994 

UVES Blue 

23.30 ±0.02 




Elliott et al. (2014) 

2012.022 

UVES Blue 

23.80 ±0.02 




Elliott et al. (2014) 


Note. — In some previous analyses, contrast ratios were not listed for specific epochs. Observations without 
listed separations correspond to simultaneous multiband photometry. 

^ Observations published without uncertainty estimates; we choose conservative values. 


ation expected from 51 Eri is 3 cm well below 

our sensitivity. 

We calculate posterior distributions for all parameters 
using emcee (Foreman-Mackey et al. 2013), an imple¬ 
mentation of the afHne-invariant Markov Chain Monte 
Carlo ensemble sampler of Goodman & Weare (2010). 
After performing a local optimization to determine a 
maximum-likelihood fit, we move 3000 walkers each 4000 
steps. We discard the first 2000 steps of each walker 
as burn-in, and use the test of Geweke (1992) and vi¬ 
sual inspection to verify the system has converged. The 
data and allowed orbits are shown in Figure 1. Summary 
statistics for the orbital parameters are given in Table 2. 
We note the fitted systemic RV of 20.76 ± 0.18 km s“^ 
is consistent with the measured RV for 51 Eri, 21.0 ± 1.2 
km s“^ (Bobylev 2006) and the UVW velocities are con¬ 
sistent with Mamajek & Bell (2014). Our samples are 
available online. 

We estimate bolometric luminosities for both stars 
by integrating the CFHIST2011_2015 model spectra of 
Baraffe et al. (2015). We use the 3700 and 3500 K mod¬ 
els with log 5 = 4.5 (cgs) as spectral templates, scaling 
them until they match the observed combined and dif¬ 
ferential magnitudes in each available bandpass. We add 
in quadrature 0.10 mag of uncertainty in our visible-light 

http://www.astro.caltech.edu/~btm/research/gj3305.html 


TABLE 2 

Parameters for G.J 3305 AB 


Parameter 

Median 


Uncertainty 

(lo-) 

y^coscj 

0.160 

± 

0.019 

-Resina; 

-0.406 

± 

0.015 

Eccentricity 

0.19 

± 

0.02 

Argument of Periastron oj [deg] 

-69 

± 

3 

Time of Periastron [Year] 

2007.14 

± 

0.16 

Orbital Period [Year] 

29.03 

± 

0.50 

GJ3305A Mass [Mq] 

0.67 

± 

0.05 

GJ3305B Mass [M©] 

0.44 

± 

0.05 

Total System Mass [Mq] 

1.11 

± 

0.04 

Mass Ratio Mq/Mj^ 

0.65 

± 

0.10 

Orbital Inclination, i [deg] 

92.1 

± 

0.2 

Orbital Semimajor Axis, a [AU] 

9.78 

± 

0.14 

Long, of Ascending Node, Q [deg] 

18.8 

± 

0.2 

Systemic RV Velocity, 7 [km s“^] 

20.76 

± 

0.18 

RV semiamplitude Ka [km s“^] 

4.01 

± 

0.38 

U [km s“^] 

-13.76 

± 

0.24 

V [km s“^) 

-16.40 

± 

0.40 

W [km s“l[ 

-9.71 

± 

0.36 

GJ3305A Luminosity [Lq] 

0.112 

± 

0.007 

GJ3305B Luminosity [Lq] 

0.043 

± 

0.005 


magnitudes and 0.05 mag in the near-IR to account for 
stellar variability. 
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Fig. 1.— (Top Left) Astrometry for GJ3305AB. Data points correspond to the observations listed in Table 1. Blue lines correspond 
to random draws from the posterior distributions of orbital elements. The red, dashed line corresponds to the maximum likelihood orbit. 
(Top Right) RV data for GJ 3305 A from the literature. The published uncertainties are in black; in gray are the best-fitting uncertainties, 
incorporating an RV jitter model. The red, dashed line corresponds to the maximum likelihood orbit. The blue shaded regions correspond 
to the 1—, 2—, and 3 (t uncertainties in the RV of GJ3305A. (Bottom Left) Measured separations for GJ3305AB and residuals from the 
maximum likelihood model. Each feature on the plot retains its meaning from the previous subplot. (Bottom Right) Measured position 
angles for GJ 3305 AB and residuals from the maximum likelihood model. 


4. GOMPARISON WITH BHAG15 EVOLUTIONARY 
MODELS 

Given the known distance to the system from Hippar- 
cos we can test if theoretical stellar evolution models 
accurately predict the inferred stellar masses and age of 
the j3 Pic moving group. Combined-light photometry 
spanning from B (0.4 /im) to Ks (2.3 /tm) was mea¬ 
sured by the APASS, 2MASS, and WISE surveys (Ta¬ 
ble 3). We add an uncertainty of 0.03 mag in quadra¬ 
ture to the listed APASS uncertainties due to the pres¬ 
ence of systematics in APASS DR9 at that level (Henden 
et al. 2012). We also have obtained one epoch of differ¬ 
ential photometry in two visible-light bandpasses with 
DSSI and two near-IR bandpasses {H and Bry) with 
Keck/NIRC2. 

We compare the observed brightness of GJ 3305 AB to 
that predicted by the BHAC15 models of Baraffe et al. 
(2015) for two stars of masses consistent with those in¬ 
ferred during our analysis as a function of age. We 
find models of 25 Myr old stars accurately predict the 
combined-light near-IR flux for these stars, although the 


TABLE 3 

Photometry for GJ 3305 AB 


Bandpass 

Source 

Magnitude 

Uncertainty 

Combined 

B 

APASS DR9 

11.94 

0.03 

V 

APASS DR9 

10.56 

0.05 

9' 

APASS DR9 

11.27 

0.03 

r' 

APASS DR9 

10.03 

0.07 

J 

2MASS 

7.30 

0.02 

H 

2MASS 

6.64 

0.05 

K 

2MASS 

6.41 

0.02 

W1 

WISE 

6.34 

0.03 

W2 

WISE 

6.21 

0.02 

W3 

WISE 

6.16 

0.02 

W4 

WISE 

6.00 

0.04 

Resolved 

A692 

DSSI 

1.89 

0.04 

A880 

DSSI 

1.17 

0.03 

AH2 

Keck/NIRC2 

1.00 

0.02 

ABry 

Keck/NIRC2 

0.92 

0.01 

AH 

Keck/NIRC2 

1.00 

0.02 


models predict brighter V magnitudes than those ob- 
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served (Figure 2). However, star B is brighter than these 
same models predict: a 25 Myr old GJ3305B would be 
significantly brighter than what is observed. Assuming 
the stars are coeval, the models then predict a mass for 
GJ 3305 B that is 20% lower than the observed mass. 

We create a simulated spectral energy distribution for 
each star, given the measured masses and the average 
age of P Pic as measured from higher-mass stars. We in¬ 
terpolate absolute magnitudes predicted by the updated 
BHAG15 models of Baraffe et al. (2015) along isochrones 
and isomass contours to predict apparent magnitudes for 
these stars in each bandpass. We find that the total re¬ 
ceived flux is lower than predicted by the BHAG15 mod¬ 
els in each bandpass. While the flux for GJ 3305 A is con¬ 
sistent with the model predictions, GJ 3305 B is fainter 
than predicted. 

Given the observed masses, we then vary the age of 
the system, assuming both stars are coeval, to determine 
which system age would be predicted by these models 
given the observed combined and differential magnitudes. 
We apply a flat prior on the age of the system, finding the 
BHAG15 models predict an age of 37±9 Myr, consistent 
with the overall age of the moving group (24 ± 3 Myr, 
Bell et al. 2015). As the system is unambiguously young, 
we can also confirm 51Erib as a planetary mass object. 

5. DISCUSSION 

We have measured the masses and orbits of 
GJ3305AB, finding both to be consistent with the 
BHAG15 models at the 1.5cr level. In the future 
GJ 3305 AB and the gravitationally bound 51 Eri Ab will 
be able to act as an isochronal test as a coeval, co- 
metallicity quadruple system spanning stellar to plan¬ 
etary mass regimes. 

The derived period of GJ 3305 (29.03 ± 0.50 year) is 
longer than the 21 year found by Delorme et al. (2012). 
The authors of that paper did not have sufficient data to 
fit all orbital parameters, so they fixed the total system 
mass to 1.3 Mq. Given our lower mass measurement, 
it is not surprising that our measured orbital period is 
longer. 


5.1. Current Limitations 

It is possible that an unseen very low-mass star or 
brown dwarf orbiting G J 3305 B could cause us to overes¬ 
timate its mass, causing the observed 20% discrepancy. 
For the system to be stable over 20 Myr, such a com¬ 
panion would have to be in a close (P < 50 day) orbit. 
The companion would then have to be in a nearly face-on 
(i < 10°) orbit to evade RV detection. Such companions 
could be found through continued astrometric monitor¬ 
ing of GJ 3305. Such a companion would not affect our 
astrometry due to its small separation from GJ 3305 B 
and would likely not affect our photometry due to its 
low luminosity relative to the other stars in the system. 

Most PMS M dwarfs have distance measurements to a 
precision no better than 5%, meaning the total mass can¬ 
not be measured to better than 15% (e.g. Shkolnik et al. 
2012). The uncertainty in the mass of GJ 3305 AB is only 
4%: the dominant source of uncertainty in this value is 
the 1% Hipparcos parallax to 51 Eri, making this sys¬ 
tem an ideal low-mass benchmark. With a Gaia parallax 
forthcoming in the next few years, parallaxes for low- 
mass PMS stars will be improved substantially. Long¬ 


term astrometric and RV monitoring of wide M dwarfs 
is essential as parallaxes are obtained over the next few 
years. 

The uncertainty in the individual mass of each star is 
dominated by the uncertainty in the Doppler semiampli¬ 
tude. While additional astrometric observations will not 
significantly improve the measured physical properties 
of GJ 3305, additional RV observations will be impor¬ 
tant. RV observations behind AO would be especially 
beneficial, as the RV from each star could be measured 
separately, instead of a flux-weighted RV centroid. 


5.2. Dynamical Effects on 51 Eri b 

GJ 3305 AB and 51 Eri Ab exist in a dynamical config¬ 
uration that may be susceptible to Kozai-Lidov oscilla¬ 
tions (Kozai 1962; Lidov 1962), as suggested by Mac¬ 
intosh et al. (2015). In this scenario, the planet-star 
binary (51 Eri Ab) interacts secularly with GJ 3305 AB, 
leading to oscillations in inclination and eccentricity of 
the planet-star sub-system. The timescale for such an 
interaction is 


r 


planet 


Mpert 


/ 0,pert 
\ ^planet 


(1 - 


( 2 ) 


where Ppianet is the orbital period of a planet with a 
semimajor axis of Opianet about a host of mass M*, 
Mpert is the mass of a distant perturber, and Opert 
and Cpert are the semimajor axis and eccentricity of the 
perturber/planet-star “binary” orbit (see e.g. Holman 
et al. 1997). 

Although we have limited information about this sys¬ 
tem, we can estimate the timescale for Kozai-Lidov cycles 
should the mutual inclination of the 51 Eri Ab system and 
(51EriAb)-(GJ3305 AB) system satisfy 140° ^ im ^ 
40°. Taking the instantaneous sky-projected separations 
as a proxy for the semimajor axes and inferred masses 
of M* = 1.75 Mq (Simon & Schaefer 2011) and Mpert = 
1.1 Mq yields a timescale of t ^ 2 x 10® yr (1 — Cpgj.^)®/^. 
Therefore, unless the eccentricity of GJ 3305 about the 
51 Eri subsystem satisfies Cpert ^ 0.9, the timescale for 
Kozai-Lidov oscillations is longer than the age of the sys¬ 
tem, so we do not expect the Kozai-Lidov mechanism to 
have had time to induce a large eccentricity or spin-orbit 
misalignment within the 51 Eri sub-system. If future ob¬ 
servations indicate non-zero spin-orbit misalignment or 
a high eccentricity for the orbit of 51 Eri b, a primordial 
origin unrelated to the distant perturbers would be sug¬ 
gested. 
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Fig. 2.— (Top) (Left) Combined-light, unresolved and (Right) differential, resolved photometry for GJ3305AB (black) compared to 
predictions (red) of the BHAC15 models as a function of age given the observed masses and parallax. The data are consistent with an 
age larger than 25 Myr. Plotted bars along the abscissa correspond to the width of each filter and are meant to guide the eye: they do 
not represent an uncertainty. (Middle left) SED for the system, assuming a 24 ± 3 Myr age and the observed masses. Combined-light 
photometry is in black and resolved photometry in purple. While the model accurately reproduces the observed flux from GJ 3305 A, it 
overpredicts the received flux from GJ 3305 B. (Middle right) Joint posterior probability distributions on the masses of the two stars, (black) 
inferred from the astrometry and RV data and (red) predicted by the BHAG15 models given the observed combined-light and differential 
photometry assuming an age of 24 di 3 Myr. Gontours correspond to the 1-, 2-, and S-a confidence regions. The BHAC15 models predict 
a mass for GJ3305B consistent with the mass inferred from the data, but underpredicts the mass of GJ3305A by 20%. (Bottom left) 
GMD showing the absolute H magnitudes and H — K colors of GJ 3305 AB compared to theoretical models. The models provide a more 
accurate fit for GJ3305A than GJ3305B. (Bottom right) Posterior probability distribution on the age of the GJ 3305 system, calculated 
by marginalizing the joint mass-age posterior over all allowed masses, assuming both stars are the same age. The BHAC15 models predict 
an age of 37 di 9 Myr; the dashed line represents the Bell et al. (2015) age of the /3 Pictoris system. 
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